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We show that the spectrum of an ensemble of two-level systems can be broadened through 're- 
setting' discrete fluctuations, in contrast to the well-known motional-narrowing effect. We establish 
that the condition for the onset of motional broadening is that the ensemble frequency distribution 
has heavy tails with a diverging first moment. We find that the asymptotic motional-broadened 
lineshape is a Lorentzian, and derive an expression for its width. We explain why motional broaden- 
ing persists up to some fluctuation rate, even when there is a physical upper cutoff to the frequency 
distribution. 
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I. INTRODUCTION 

Ensembles consisting of many two-level system (TLS) 
are of interest in many physical disciplines. In precision 
metrology, the energy difference between such two lev- 
els in a cesium atom is used to define the second. In 
quantum computation, the TLS is called a qubit and it 
replaces the classical bit with an added ability to store 
any superposition of the two underlying logical values. 
An ensemble of such TLS can serve as a storage medium 
for quantum information lj. A quantum memory is a 
necessary building block in quantum networks since it 
relieves the requirement for a series of successful quan- 
tum operations and therefore enable scalability [H [3] . In 
order to extend the coherence time of such a memory it 
is first essential to understand the effect of its coupling 
to the environment. 

The difference in the energy levels, or the transition fre- 
quency between the two, is in general not identically the 
same for all the TLS in the ensemble, and its distribution 
has a certain breadth. For instance, for a cold atom en- 
semble held together by a confining potential, the differ- 
ence between the energy levels of the internal states, may 
depend on the location of the atom within the trap, and 
may also be affected by inter-particle interactions. Al- 
though the "bare spectrum" due to the inhomogencitics 
within the ensemble is usually time- independent, the fre- 
quencies of the individual TLS are in general not con- 
stant, and undergo time-dependent fluctuations. 

Such fluctuations are usually thought to induce nar- 
rowing of the power spectrum in observations of the de- 
cay rate of the ensemble coherence, a phenomenon named 
motional narrowing. Historically, motional narrowing 
was first observed in NMR, where the spectra of liquid 
materials were found to be significantly narrower than 
those of solids due to thermal motion of the nuclei [I]. 
Spectral narrowing due to fluctuations was later encoun- 
tered in many other fields including in hot atomic vapors 
(Dicke narrowing) [5], semiconductor microcavities [6], 
quantum dots [7] and trapped cold atomic ensembles [5] . 



In this paper we show that fluctuations can have the 
reverse effect and lead to broadening of the spectrum 
(motional broadening). In terms of quantum informa- 
tion, motional broadening manifests itself as a shorten- 
ing of the coherence time as the the fluctuation rate in- 
creases. We prove that the condition for the emergence 
of motional broadening is that the ensemble frequency 
distribution will have heavy tails with a diverging mean. 
An example for this effect was first pointed out in Ref. 
[5]. We also show that for both motional narrowing and 
broadening the asymptotic decay of the coherence is ex- 
ponential, and derive an expression for the decay rate. 
Since in practice heavy tails of the frequency distribu- 
tion can be sustained only up to some point, we study 
scenarios with cutoffs and show that motional broadening 
persists up to some fluctuation rate. The motional broad- 
ening phenomenon should be relevant to many fields in 
which heavy-tail distributions are encountered, including 
turbulence [TU], diffusion [TT] and laser-cooling [12] . 

This paper is organized as follows. In section [TTJ we de- 
scribe the model and the quantity which is measured in 
the experiment, i.e. the coherence function. To show that 
the effect of the increase in the fluctuation rate (in time) 
depends on the tail of the probability distribution, we 
shall first consider it within the context of the so-called 
'stable laws', whose defining property is naturally suited 
for our purpose (section III). The rule is then shown to 



apply also beyond those special cases by first demonstrat- 
ing it on a family of probability distributions with a tun- 
able parameter (the 'student's t-distribution', section |v|. 
The main result is then generalized to distributions which 
converge to stable distributions (section VI). 



Further discussion is presented in section | VII| we note 
that motional broadening is analogous to the anti-Zeno 
effect, and draw an analogy between the spectroscopic 
problem we have considered and diffusion in real space. 
We use this analogy to propose an experiment to observe 
the phenomenon. We continue with the discussion of the 
effect of a cutoff in the distribution, and the asymptotic 
shape and width of the broadened spectrum. Our con- 
clusions are given in section |VIII| 
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II. THE MODEL 

We consider an ensemble of two- level systems (TLS) 
with internal states designated by |1) and |2), which are 
eigenstates of the Hamiltonian. Their energies are of the 
form 



E h2 = Eo±h(u + 6(t))/2 



(1) 



with Tiljo the 'base line' energy difference, and Ti5(t) the 
detuning from this base line. S(t) is a random correction 
whose value fluctuates within the ensemble and which 
for individual TLS also changes discontinuously in time. 
We assume a steady-state situation in which the distribu- 
tion of the 'detuning' term 5(t) is time-independent, and 
denote by Po(S)dS its distribution over the ensemble. 

In a Ramsey experiment - like procedure, e.g., as de- 
scribed in 8 , a tt/2 pulse may be used to initiate the 
ensemble's individual TLS in the coherent superposition 
states \tp) — l/\/2(|l) + |2)). In a situation where the 
internal eigenstates can be treated as stationary, but the 
energy differences between them fluctuate in time the in- 
dividual TLS states' density matrices evolve into 

\HT)) (V(t)l = ^[|l)(l| + |2) (2|] 

+ l e -i[w T+4>(T)]/2 ijv , 2 | + 1 e i[uj T+4>(T)}/2 ^ ^| _ 

(2) 

where the fluctuating part of the accumulated phase dif- 
ference between the two internal states is 



4>{T) = / 5{t)dt. 
Jo 



(3) 



The phase difference in the off-diagonal terms can be 
determined using a second tt/2 pulse, applied at time 
t. The measured quantity is described by the coherence 
function: 



R{T) = (e 



i$(T)\ 



(4) 



where (...) denotes the ensemble average [IB] . 

Our discussion concerns the effect on the ensemble co- 
herence function of the rate at which the individual de- 
tuning terms S(t) are refreshed. In a model which may 
capture the effect of 'hard collisions' these terms change 
discontinuously at random times, {tj}, at which the value 
of 5(t) is reset with the stationary distribution P (5)d5. 

The coherence function starts at i?(0) = 1 and decays 
to on a time scale which is referred to as the coherence 
time. If the absolute value is omitted from its definition, 
oscillations may occur. The spectrum function, S(lj), is 
the absolute value squared of the coherence function's 
Fourier transform. In general, the width of the spectral 
distribution varies in the opposite way to the length of 
the coherence time, e.g, spectral narrowing corresponds 
to the extension of the coherence time. 



As long as the fluctuations in time are negligible, the 
detuning of each TLS occurs at a constant rate and the 
ensemble coherence is given by 



Ro(T) = 



P (6)e lST d6 



(5) 



It was already noted that under certain conditions the 
coherence time is increased due to fluctuations (motional 
narrowing). The main result presented here is that the 
coherence time can also be reduced by the increase of the 
fluctuation rate. The latter happens when the detuning 
distribution Po(S) has heavy tails, with a diverging first 
moment (to avoid confusion, we shall henceforth reserve 
the term fluctuations for fluctuations in time of individ- 
ual TLS, not to the statistical fluctuations which occur 
within the ensemble). 



III. MOTIONAL BROADENING FOR STABLE 
DISTRIBUTIONS 

It is instructive to consider the effect of fluctuations 
when the distribution of the detuning S is given by one of 
the so-called 'stable laws' with a characteristic exponent 
< a < 2 [16]. Some well-known examples of stable 
distributions are Gaussian (a = 2), Cauchy (a — 1) and 
Levy (a = 1/2) distributions. For each distribution in 
this class, the weighted sum of independent identically 
distributed variables produces a variable with a scaled 
version of the same distribution. More explicitly: for 
any two real numbers q, s > 0, and a pair of independent 
variables 5\ , 82 of such distribution, the weighted sum 
(qSi + 5^2) has the same distribution as (q a + s") 1 '™ S. 
The characteristic function of any a-stable distribution 
satisfies [16] : 



\Mt)\ 



-itS\ 



(6) 



with some c a > 

Assuming P is a-stable, the accumulated phase due 
to the discrete fluctuation events ('collisions') at which 
6(t) is reset can be written as 



l/a 



4>{T) 



j=i 



\i =1 



5T 



(7) 



with Atj denoting the periods between collisions, Tj 



D 



Atj/T, and = standing for the equivalence of the cor- 
responding distributions. By (|6| the coherence without 
collisions is given by Rq(T) = e~ CaT and with collisions 
it is given by 

R(T) = e - c ^=i At ?) = ,r ( T )£?=i t " . (8) 
For any series of collisions: X)j=i T j = 1) ano - hence: 



, . < (w)"- 1 a > 1 
" 1 (9) 

> (w)"" 1 a < 1 
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with r max = max{r,'}. 

Combining Q , ^ and the fact that T max < 1 we find 
that the transition from motional narrowing to broaden- 
ing occurs at a = 1: 

R(T) I > R °( T ^ ' a> 1 ( motional narrowing) 
1 < Ro(T) , a < 1 (motional broadening). 

It may be noted that: 

1. The above holds regardless of the distribution, or 
the values, of the collision times 

2. If the time between collisions is constant, Atj = 
t co u, then the coherence function decays exponen- 
tially in T: 

R(T) = e-MW- 15 ^ (11) 

That is: the coherence time is given by 
Ca(icoiz)^" -1 ^ which is in line with the observation 
made in (10). A similar conclusion holds if Atj are 
random but of common order t co u. 



IV. THE SOLVABLE CASE OF POISSON 
FLUCTUATIONS 

An explicit expression for R(t) can be derived for dis- 
crete fluctuations (e.g., collisions) which occur with a 
Poisson distribution in time, at rate T. After a random- 
izing event, the TLS acquires a new detuning rate, inde- 
pendently drawn with the probability distribution Pq(5). 
The coherence function can be calculated exactly in this 
model through the following expression for its Laplace 
transform, R(s) = / °° e~ sT R(T)dT, [HEI]: 



r = u.j 






r/<s = o 


r/6 = io 
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FIG. 1. The two upper graphs show the sp ectr um for a Stu- 
dent's t-distribution of the detunings [see ( 13 1] . On the left 
the spectrum is plotted for r — 0.5 for which the first moment 
of the distribution diverges whereas on the right the spectrum 
is plotted for r = 1.5 for which only the second moment di- 
verges but the first moment exists. For both spectrums the 
spectrum is plotted without fluctuations (solid line) and with 
fluctuations at T — KWo (dashed line). In the lower graph we 
plot the full width at half the maximum (FWHM) normalized 
to the FWHM without fluctuations, as a function of T. 



with the normalization factor 



Y,r.V) = r(^J/r(£)W^ , (14) 



R(s) 



Ro(s 



i -ri? (s + r) 



This relation is valid for all Po(S). It was shown ex- 
perimentally that ( 12 ) correctly describes the effect of 
cold-collisions in trapped atomic ensembles |15j . Equa- 
tion ( 12 ) is also useful for model calculations, such as the 



one presented next. 



V. AN ILLUSTRATIVE EXAMPLE 

To illustrate the transition from motional narrowing 
to broadening beyond the above class of special distri- 
butions, we consider an ensemble with a Student's t- 
distribution, with a tunable parameter r: 



P (5)=N(r,6 ) 



1/5 

r U 



(12) where T(z) is the gamma function. For r — > oo the dis- 
tribution is approaching a Gaussian with a standard de- 
viation 6q, and for r = 1 it is identical to the Cauchy 
distribution (Lorentzian). The first (second) moment 
of the distribution diverges for r < 1 (r < 2). Using 
(12 1, we calculate the spectrum and plot it in Fig. [ljfor 
r = 0.5 and r = 1.5, with and without fluctuations. For 
r = 1.5 we observe that the spectrum becomes narrower 
in the presence of fluctuations, which demonstrate that 
motional narrowing persists even when the second mo- 
ment diverges. On the other hand, for r = 0.5 the fluctu- 
ations broaden the spectrum. In Fig. [T]we plot as a func- 
tion of r the normalized spectral width, which is defined 
to be the full width at half the maximum (FWHM) di- 
vided by the FWHM for T = 0. The figure clearly shows 
the narrowing (for r — 1.5) or broadening (for r = 0.5) 
effects as the fluctuation rate increases. Curiously, for 
the Cauchy distribution, corresponding to r — 1, there is 
^3) no r dependency This fact is in line with both (10) and 



(12) 
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VI. GENERALIZATION 

The above observations can be extended further to dis- 
tributions which are by themselves not stable, but are in 
the domain of attraction of a stable distribution (or 'law') 
S a . This notion means that a sum of variables drawn 
from the distribution, up to a normalizing factor, con- 
verges in distribution to S a , as the number of summands 
increases. 

A distribution belongs to the domain of attraction of 
an a-stable law if its cumulative distribution function, 
F(x), scales as F(x) ~ |x| Q /i(|a;|) as x — > — oo, and 1 — 
F(x) ~ x a h(x) as x — > oo, where h(x) is a slowly varying 
function at infinity [17] . The domain of attraction of 
the Gaussian distribution contain all distributions with 
a finite variance. 

If Pq belongs to the domain of attraction of an a- sta- 
ble distribution, then its characteristic function is of the 
form \ip4,{t)\ = e^l'l ^*), with h(t) = e o(t) a slowly vary- 
ing function as t — > [17]. The coherence without colli- 

and with col- 



sions is given by Rq(T) = exp —cT a h(T) 



lisions it is given by R(T) — exp 
Extending ([9| one may see that if lim ra 
fixed T, than for a > 1: J2"=i Atfh{Atj 
and for a < 1: £™ Atfh(Atj) > C 2 T a T m 2l, with some 



J = J- 3 
oo Tmax — 0? & 
\ ^i-i 'mux' 



rjj , Cj > 0. This extends the validity of ( 10 1 for the limit 
of many randomizing events (high collision rate) when 
the detuning have a distribution in the domain of attrac- 
tion of an a-stable law. 



VII. DISCUSSION 

The above can be summarized by saying that whether 
the coherence of a TLS ensemble decays faster or slower 
due to 'resetting' discrete fluctuations depends on the the 
tails of the detunings distribution. Motional broadening 
emerges for heavy-tailed distributions with a < 1 which 
corresponds to a diverging first moment. This explains 
the results of Fig. [T] since the Student's t-distribution be- 
longs to the domain of attraction of an a-stable distribu- 
tion with a = r. Furthermore, the criterion for motional 



broadening given in ( 10 ) coincides with the divergence of 



the detuning distribution's first moment. 

The relation to the anti-Zeno effect. — To get an intu- 
ition for the above results we write the coherence at a 
time T, for a given series of randomization events, as 



n „oo 

R(T)=l[ ^P (^K AM 
i=i J -°° 



RoiAh) ■ R {At 2 ) ■ ■ ■ R (At n ) 



(15) 



In addition, the derivative 9t-Ro|t=o+ is for a > 1 and 
— oo for a < 1. The combination of these two properties 
explains why depending on whether a is larger or smaller 
than 1, the coherence after a resetting collision lie above 



r = 1.5 

motional narrowing 



r = 0.5 
motional 
broadening 



Time [So 1 ] 



FIG. 2. The coherence without fluctuations (dotted lines) and 
with fluctuations separated by Atj = 0.25 ■ 6q (solid lines). 
The detuning distribution is assumed to be the Student's t- 
distribution with r = 0.5 and 1.5. Since the coherence is given 
by R(T) = iio(Ati) ■ Ro(At 2 ) ■ ■ ■ R (At„), a Zeno/anti-Zeno 
like behavior explains the transition from motional narrowing 
to broadening for a diverging first moment of Po, at which 
point 8tRo\t=o changes from to — oo. 



or below the original curve of Rq(T), as depicted in Fig. [2] 
for the simple case of equal times between such collisions. 
In this respect, motional narrowing is analogous to the 
Zcno effect [IS], and motional broadening to the anti- 
Zcno effect. 

Asymptotic lineshape. — We rewrite ^ using the typi- 
cal time between collisions Atj ~ r , and the inhomo- 

geneous decay rate 70 = c« , and obtain in the limit 
of many collisions FT > 1: R(T) w e -^^- a T _ The 
asymptotic behavior of the coherence decays exponen- 
tially with time, and the decay rate is given by 



7 = 7o r 1 



(16) 



This equation for a — 2 was verified experimentally with 
optically-trapped cold atomic ensemble [S]. Since Q is 
true only for a stable distribution, we test the validity 
of ( 16 ) for distributions in the domain of attraction of 
an a-stable distribution in numerical simulations. The 
results of these simulations for a Student's t-distribution 
are plotted in Fig. [3] and show that as the fluctuation 
rate increases the decay of coherence indeed becomes ex- 
ponential. From the numerical curves we extract the de- 
cay rate and plot it in the inset of Fig. [3] for various 
values of the distribution parameter r and T. The func- 
tional form of the decay rates confirms the prediction of 



( 16 1 with the correct a. 



The effect of a cutoff. — In real physical situations the 
detuning distribution can not have a diverging first mo- 
ment, and the heavy tail scaling can be sustained up to 
some cutoff 5 C . For Po(5) with a characteristic exponent 
a < 1, the order of magnitude of the sum J2j=i ^tjSj 
is the same as of max{^<5j}. The effect of the cutoff is 
therefore negligible as long as Prob(max[{(5j}] > S c ) <C 1 
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FIG. 3. Numerical simulation of the logarithm of the co- 
herence using an ensemble of 10000 particles following the 
Poisson discrete fluctuations model with three different rates 
r and detunings following a Student's t-distribution with 
r — 0.5. The dotted (black) lines are linear fits, validating 
the prediction of (161 that as F increases the decay becomes 
exponential. We extract the decay rate, 7, by fitting the sim- 
ulated coherence for T > lOr" 1 to an exponentially decaying 
function Ae fT . The inset shows 7 as a function of T for dif- 
ferent values of the distribution parameter r. The dotted lines 
are fits to the function 7 — aF 1 '" + b (which has the func- 
tional form derived in (16 1. The extracted exponents agree 



well with the expected values a = r. 
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FIG. 4. The effect of a cutoff in the detuning distribution on 
the motional broadening phenomenon. The figure shows nu- 
merical simulation of the normalized coherence R(T)/Rq(T) 
at time T — O.HSg 1 versus the fluctuations rate F. The de- 
tuning distribution is taken to be the Student's t-distribution 
[see < \13\ ] with r = 0.5, and the fluctuations are Poissonian. 
Each curve is calculated for Po truncated at a cutoff detuning 
S c . As expected in motional broadenings, the coherence is de- 
creasing for small V, and only for higher values the behavior 
is changing to motional narrowing. 

(we assume the average (At,-) = < 00). This 

probability depends on the number of collisions, which 



is roughly given by TT. An estimate of this probabil- 
ity yields that the effect of the cutoff is insignificant for 
TT <C {S c /So) a , where S is the typical scale of the detun- 
ing distribution [see (13)]. This means that for a given 
observation time T, motional broadening persists up to 
fluctuation rate on the order of (<5 c /(5o)"T _1 . This qual- 
itative picture is demonstrated in numerical simulations 
plotted in Fig. [4] Motional broadening prevails for small 
r, later changing to motional narrowing once the cutoff 
is sampled and discovered. 

The analogy to diffusion in real space. — The TLS en- 
semble coherence problem can be mapped to that of par- 
ticles performing diffusion in real-space, where the detun- 
ing rate and accumulated phase are mapped to velocity 
and position, respectively. In this analogy, the diffusion 
problem assumes an ensemble of particles with a steady- 
state distribution of velocities, starting all from the same 
point in space. In the absence of collisions, the particles 
are ballistically expanding and the width of their position 
distribution grows linearly with time. With collisions, 
our criterion yields that the width of the particles' spatial 
distribution grows faster than ballistic (super-ballistic) 
for heavy-tailed velocity distribution. Super-ballistic dif- 
fusion is known to exist in turbulent flow [19] . In the 
context of spatial diffusion, a particularly interesting im- 
plementation of the model considered in this paper can 
be achieved. It was shown that the steady-state velocity 
distribution of atoms in a polarization lattice follows a 
power law with an exponent that depends on the lattice 
depth [20] . As a result, in such a system the diffusion 
becomes anomalous 21. Going to a low enough lattice 
depth, it may be possible to observe motional broaden- 
ing, namely diffusion whose scaling with time is faster 
than ballistic. 



VIII. CONCLUSIONS 

In conclusion, inhomogeneous broadened spectrum of 
an ensemble of two-level systems is modified when the 
energy of each system is fluctuating. We have shown 
that although most commonly the fluctuations lead to 
narrowing of the spectrum, under certain conditions the 
opposite can occur. This crossover is determined by the 
tail behavior of the inhomogeneous energy distribution, 
with motional broadening arising for distributions of di- 
vergent first moment. Finally we have drawn the analogy 
to the anti-Zeno effect and to faster-than-ballistic spatial 
diffusion. 
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